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■^j- ' Abstract. Angular power spectra and structure functions of the Stokes parameters Q and U and polarized intensity 

, P are derived from three sets of radio polarimetric observations. Two of the observed fields have been studied at 

\Q ■ multiple frequencies, allowing determination of power spectra and structure functions of rotation measure RM 

as well. The third field extends over a large part of the northern sky, so that the variation of the power spectra 
over Galactic latitude and longitude can be studied. The power spectra of Q and U are steeper than those of 
P, probably because a foreground Faraday screen creates extra structure in Q and U, but not in P. The extra 
structure in Q and U occurs on large scales, and therefore causes a steeper spectrum. The derived slope of the 
power spectrum of P is the multipole spectral index cep, and is consistent with earlier estimates. The multipole 
|' spectral index otp decreases with Galactic latitude (i.e. the spectrum becomes flatter), but is consistent with a 

1~ \ constant value over Galactic longitude. Power spectra of the rotation measure RM show a spectral index a am ~ 1, 

while the structure function of RA1 is approximately flat. The structure function is flatter than earlier estimates 
from polarized extragalactic sources, which could be due to the fact that extragalactic source RM probes the 
complete line of sight through the Galaxy, whereas as a result of depolarization diffuse radio polarization only 
r* ■ probes the nearby ISM. 
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1. Introduction Small-scale structure in the warm ISM and magnetic 

field can be well studied using radio polarimetric ob- 

The warm ionized gaseous component of the Galactic in- servationSj of the synch rotron background in the Milky 

terstellar medium (ISM) shows structure m density and Way (e g g rouw & Spoelstra 1976, Wieringa et al. 1993, 

velocity on scales from AU to several kiloparsecs. The Duncan et aL 1997; 1999; Uyamker et al. 1999, Landecker 

Galactic magnetic field is coupled to the motions of the et aL ^ Gaensler et aL 2 001), of pulsars (e.g. Rand & 

ionized gas and has a comparable energy density, so that Kulkarni 1989; hno & Shibata 1993, Rand & Lyne 1994, 

gas and magnetic field are in complex interaction. Detailed Han et al 19Q9) or of polarized extragalactic point sources 

knowledge of the turbulent nature of the warm ISM and (e g Simard-Normandin & Kronberg 1980, Clegg et al. 

the structure in the Galactic magnetic field is essential for 1992) ^ At ghort wave i e ngths ( A < 6 cm) , Faraday rotation 

several fundamental studies of the Galaxy, including mod- is negligible; so that the measured polarization directly 

eling of molecular clouds (e.g. Vazquez-Semadeni & Passot traces the magnetic field in the emitting region. At longer 

1999, Ostriker et al. 2001), heating of the ISM (Minter & wave i e ngths, Faraday rotation measurements give addi- 

Balser 1997), star formation (e.g. Ferriere 2001), and cos- tional information on dens it y an d magnetic field structure 

mic ray propagation (Chevalier & Fransson 1984). along the entire Une of sight Furtherm ore, depolarization 

processes define a distance beyond which polarized radia- 
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tion is significantly depolarized, which depends on wave- 
length. So high-frequency measurements probe the total 
line of sight through the Galaxy, whereas low-frequency 
polarization observations only trace the nearby part of 
the ISM. 

Specific intriguing small-scale structures and discrete 
objects have been studied in diffuse polarization observa- 
tions (Gray et al. 1998, Haverkorn et al. 2000, Uyamker 
& Landecker 2002), but the spatial structures have also 
been analyzed statistically. Simonetti et al. (1984) and 
Simonetti & Cordes (1986) studied the structure in 
Galactic rotation measure RM by comparing RMs of po- 
larized extragalactic sources, and of separate components 
of the same source. The electron density seems to exhibit 
a power law structure function and therefore a power law 
angular spectrum (see also Armstrong et al. 1995, Minter 
& Spangler 1996). Recently, statistical analysis of the dif- 
fuse Galactic polarized foreground have been pursued in 
the form of angular power spectrum studies (Tucci et al. 
2000, 2002, Baccigalupi et al. 2001, Giardino et al. 2002, 
Bruscoli et al. 2002), with the objective of estimating the 
importance of the Galactic ISM as a foreground contam- 
inator for CMBR polarization observations (e.g. Seljak 
1997, Prunct et al. 2000). 

In this paper, we study the statistical properties of the 
warm ISM and Galactic magnetic field by means of power 
spectra of Stokes parameters Q and U , polarized intensity 
P and rotation measure RM. We also derive the structure 
function of RM to allow a comparison with earlier studies 
of RM structure functions from polarized extragalactic 
sources. Furthermore, by careful selection of reliable RM 
determinations (i.e. those with low \ 2 of the linear fit to 
4>(X 2 )) in the calculation of structure functions, we can 
obtain an estimate of how much the structure functions 
(and power spectra) are influenced by low-quality RMs. 

We use data from three regions, all at positive Galactic 
latitudes, in which we observed the diffuse polarized emis- 
sion at frequencies around 350 MHz with the Westerbork 
Synthesis Radio Telescope (WSRT). The first two regions, 
each about 50 square degrees in size, in the constella- 
tions Auriga (Haverkorn et al. 2003a) and Horologium 
(Haverkorn et al. 2003b), are observed at multiple frequen- 
cies at a resolution of ~ 5'. The third region is a part of 
the Westerbork Northern Sky Survey (WENSS, Rengelink 
et al. 1997), a high-resolution radio survey at 327 MHz of 
the northern hemisphere (Schnitzeler et al., in prep). For 
those parts of the WENSS survey that were observed at 
night, polarization data are usable. Here we discuss po- 
larization data from the region with 140° ^ I Ss 170° and 
0° <S b <, 30°. 

In Sect. |21 we describe the three sets of data that we 
analyze in this paper. In Sect. angular power spectrum 
analysis is introduced, power spectra of the data are pre- 
sented and discussed and literature on the angular power 
spectra is briefly summarized. Section 21 gives structure 
functions for RM in the two multi-frequency measure- 
ments. In Sect. the results are discussed, and finally in 
Sect. El some conclusions are stated. 



2. The observations 

2.1. Multi-frequency WSRT observations 

We carried out low frequency radio polarimetry with the 
Westerbork Synthesis Radio Telescope (WSRT), in two 
regions of the sky at positive latitudes, in the constel- 
lations of Auriga and Horologium. Data were obtained 
simultaneously at 8 frequencies around 350 MHz, each 
with a bandwidth of 5 MHz. Due to radio interference 
and hardware problems, only data in 5 frequency bands 
could be used, viz. those centered at the frequencies 341, 
349, 355, 360, and 375 MHz. The multi-frequency data 
allow the study of the frequency dependence of the polar- 
ization structure, and the determination of rotation mea- 
sure RM. We use the technique of mosaicking (i.e. the 
telescopes cycle through a number of adjacent fields on 
the sky during a 12hr observation) to obtain a field of 
view that is larger than the primary beam. Mosaicking 
also suppresses instrumental polarization to below 1% (see 
Haverkorn 2002). Maps of the Stokes parameters /, Q, and 
U were derived from the observed visibilities. RMs were 
computed straightforwardly from the linear relation be- 
tween polarization angle <f> and A 2 . Because the observed 
RM values are small (|i?M| Ss 10 rad m~ 2 ), there is no 
n 180° ambiguity in cf>, and RMs can be computed with 
\4>(Xi) — 4>(^j)\ < 90° for adjacent wavelengths Xi and Xj. 
We define a determination of RM in a particular beam 
reliable if (1) the reduced x 2 of the linear (/>(A 2 )-relation 
Xred < 2, and (2) the polarized intensity averaged over 
wavelength (P) > 20 mJy/beam (i.e. ~ 4<r). 

The maximum baseline of the observations was 2700m, 
yielding a resolution of 1', but smoothing of the Stokes Q 
and U data (using a Gaussian taper in the (u, w)-plane) 
was applied to obtain a better signal-to-noise. The taper 
has a value of 0.25 at a baseline value around 300m, where 
the exact values were chosen so that the beam size is iden- 
tical at all 5 frequencies, viz. 5.0' x 5.0' cosec 8. 

The first region, in the constellation Auriga, is cen- 
tered on (l,b) = (161°, 16°) and is about 9° x 7° in size. 
The left panel in Fig. ^ shows the polarized intensity P 
at 349 MHz in the Auriga region at 5.0' resolution. We 
do not show the original 1' resolution map, because it is 
noise-dominated. The maximum polarized brightness tem- 
perature in the map is X^poZ ~ 13 K, and the noise is 
about 0.45 K. Rotation measures in the Auriga region are 
shown in the right panel of Fig. ^ as circles superimposed 
on the grey scale map of P at 349 MHz. The diameter of 
each circle indicates the value of the RM at that position, 
where filled circles denote positive RMs, and open circles 
negative RMs. We show only the RMs that are reliably 
determined, and only one in four independent beams. 

The region in the constellation Horologium, centered 
at (I, b) = (137°, 7°), was observed in the same way as the 
Auriga region, at the same five frequencies, and with the 
same taper applied. Fig. |3 shows the polarized intensity 
P at 349 MHz in the left panel, and RMs in the form of 
overlaid circles in the right panel. RMs are coded in the 
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Fig. 1. Grey scale representation of the polarized intensity P at 349 MHz and rotation measure RM, of the region 
in the constellation of Auriga centered at (l,b) = (161°, 16°) at a resolution of ~ 5'. The grey scale is saturated at 
P = 95 mJy/beam (white) (Tb. po i = 12 K). The right panel shows RM in superimposed circles, where the diameters 
of the circles indicate RM, and filled circles denote positive RMs. The scaling of RM on the right-hand side is in 
rad m~ 2 . 




55 50 45 40 55 50 45 40 

a 1950 (°) a 1950 (°) 



Fig. 2. P and RM for the region in the constellation of Horologium centered at (l,b) = (137°, 7°). Notation as in 
Fig.^ where the maximum P = 95 mJy/beam corresponds to Tb tPO i — 13.9 K. 



same way as in Fig.^ The maximum polarized brightness 
temperature is Tb iPO ; ps 17 K, and the noise is a little higher 
than that in the Auriga field, about 0.65 K. 

The polarized intensity is not corrected for noise bias 

Pdebias = ViQ'tbs + U"ab») ~ °* ( for P ><>)■ However, P 
generally has a S/N > 4 — 5, for which the debiasing does 
not alter the data by more than 2 - 3%. Furthermore, 



power spectra are not affected by debiasing except at zero 
frequency. 

In both fields there is no small-scale structure visible in 
the map of total intensity Stokes /, despite the ubiquitous 
structure on arcminute and degree scales in P. Because 
the large-scale (Jj 1°) component of I cannot be measured 
with the WSRT due to missing short spacings, we used the 
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Haslam ct al. (1981, 1982) radio survey at 408 MHz to es- 
timate the total intensity I to be ~ 34 K and ~ 47 K in 
Auriga and Horologium, respectively. The polarized inten- 
sity P shows structure up to ~ 10 - 15 K on scales from 
arcminutes to degrees. Due to the lack of corresponding 
structure in /, the structure in polarized intensity can- 
not be due to small-scale synchrotron emission but must 
be due to other, instrumental and depolarization, mecha- 
nisms. 

In a medium that emits synchrotron radiation and si- 
multaneously causes Faraday rotation, the polarized emis- 
sion is depolarized by so-called depth depolarization, due 
to the vector averaging of contributions from different 
parts of the line of sight. This, together with beam depo- 
larization (due to angle structure within one synthesized 
beam), produces structure in P. Furthermore, the insen- 
sitivity of the interferometer to large-scale structure can 
cause additional structure in P. However, we have shown 
that this effect cannot be very important in these obser- 
vations (Haverkorn et al. 2003c). 

2.2. Polarization data from the WENSS survey 

The Westerbork Northern Sky Survey (WENSS, 
Rengelink et al. 1997) is a low-frequency radio sur- 
vey that covers the whole sky north of 5 = 30° at 
325 MHz to a limiting flux density of approximately 
18 mJy (5cr ) and with a resolution of 54"x54" cosec 5. 

Polarization data taken during the day are greatly 
affected by solar radiation, which is detected in side- 
lobes. Furthermore, ionospheric Faraday rotation rapidly 
changes during sunrise and sunset, causing a reduction of 
the apparent polarized intensity, and adding much noise. 
However, in mosaics observed (almost) entirely during 
night time, the polarization data are of good quality. As a 
result, we could make a "supermosaic" of a large region of 
- 30° x 35°, which we will refer to as the WENSS polar- 
ization region. Details on the observations and correction 
for ionospheric Faraday rotation are given in Schnitzeler 
et al. (in prep). 

Fig. G] shows polarized intensity P in the WENSS po- 
larization region, resampled in Galactic longitude and lat- 
itude. For this analysis, a Gaussian taper with value 0.25 
at a baseline of 500m was applied, which yields a resolu- 
tion of ~ 2.5'. In the figure, P saturates at 35 mJy/beam, 
which coincides with a polarized brightness temperature 
of ~ 17 K to ~ 25 K, depending on declination. The aver- 
age P ps 2.6 mJy/beam. Note that these observations are 
taken in a single 5 MHz wide frequency band, so rotation 
measure data are not available. 

3. Angular power spectrum analysis 

3.1. Determination of the multipole spectral index 

To quantify the structure in the polarization maps, we 
calculate angular power spectra PS(£) as a function of 
multipole £. A multipole £ is a measure of angular scales 



equivalent to wave number, and is defined as £ ps 18O°/0, 
where 9 is the angular scale in degrees. The angular power 
spectrum PS of a radiation field X is the square of the 
Fourier transform of X: PS X {£) = |J^(X)| 2 , where T de- 
notes the Fourier transform. Here, the observable X can 
be either Stokes Q, Stokes U, polarized intensity P or ro- 
tation measure RM . The power spectra were computed in 
two dimensions, and averaged over azimuth in radial bins. 
The multipole spectral index a, defined as PSx(£) oc £~ a , 
is calculated from a log-log fit to the power spectrum. 
In the tapered data, multipoles with higher values of £ 
are affected by the tapering, and in the untapered data 
higher multipoles are dominated by noise. Multipoles with 
£ ^ 200 correspond to angular scales 9 J; 1°, to which the 
WSRT is not sensitive. 

The visibilities from which a map is made are 
V map (u,v) = V obs (u,v)T(u,v), where V obs are the ob- 
served visibilities, and T(u,v) is the taper function. The 
calculated intensity of the tapered data is J-(V ma p) = 
3~(V bs) * 3~{T)i where J 7 is a Fourier transform and the 
asterisk denotes convolution. The power spectrum of the 
Stokes parameter X, PS(X), is then: 

PS X (£) = \T(X)\ 2 = \T(X obs )\ 2 T 2 = PSx, obs (t)T 2 (1) 

where X is Stokes Q or U. Although polarized intensity P 
is derived from Q and U and thus not directly observed, 
correction for the taper in the same way as for Q and U 
power spectra is a good approximation. As an illustration, 
Fig. 0] shows the power spectrum of P of the tapered data 
in the Auriga region at 341 MHz (solid line). The dotted 
line is the same spectrum, but corrected for the tapering 
according to Eq. (JIJ. The power law behavior extends to 
log(f) ps 3.6. 

3.2. Power spectra from the multi-frequency WSRT 
studies 

In Fig. we show the power spectra of P in the Auriga 
and Horologium regions at 5 frequencies, both for the ta- 
pered (upper curve) and untapered data (lower curve) 
in the same plot. The amplitudes of the power spectra 
of the tapered data are lower than those of the unta- 
pered data only because the intensities are expressed in 
mJy/beam. Because the beam widths are different for the 
two datasets, this gives a difference in the magnitude of P 
in tapered and untapered data. The power spectra of RM 
in the Auriga and Horologium field are given in FigEl Only 
tapered data give reliable enough RM determinations to 
produce power spectra for them. 

Figs. and [S] show power spectra for the Stokes pa- 
rameters Q and U in the Auriga and Horologium region 
respectively, again for tapered and untapered data. The 
corresponding multipole spectral indices a, derived for a 
multipole range of 400 < t < 1500, are given in Tabled 

At small scales (large £), the power spectra of the un- 
tapered data flatten out due to the noise in the maps, 
while the low-resolution data steepen due to the taper- 
ing, as illustrated in Fig. 0] At large scales, the Q and U 
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Fig. 3. Grey scale representation of the polarized intensity of the polarization part of the WENSS survey at 325 MHz, 
where white denotes a maximum intensity of 35 mJy/beam (Tb^ po i « 17 K to 25 K, depending on declination). The 
data were smoothed with a 500m taper, and have a resolution of ~ 2.5'. 



power spectra of the tapered data show a decrease. This 
decrease could be due to the lack of large-scale structure 
(see Sect. I2.l|) . but then it is hard to explain why there is 
no such decline in P. 

The power spectra of Q and U are steeper and have a 
larger amplitude than the power spectra of P. This could 
be caused by the presence of a Faraday screen in front of 
the emitting region. A Faraday screen will rotate the po- 
larization angle, and so induce extra structure in Q and 
U, while leaving P unaltered. This results in a higher am- 
plitude of the power spectrum. As the Faraday screen con- 
sists of foreground material, its angular size is large, steep- 
ening the spectrum. This effect was also noticed by Tucci 
et al. (2002). 

The logarithmic slope of the power spectrum of po- 
larized intensity is ap sa 2.1 - 2.3 (Table [I}, which is 
slightly higher than most earlier estimates from the liter- 
ature (although those are taken at higher frequencies, see 
Sect. l3.4|> . However, note that the Auriga and Horologium 
regions were selected for their conspicuous structure in P, 
so we expect these regions to show more structure on large 
(degree) scales than the "average" ISM, and thus exhibit 



a steeper spectrum. The power spectra of Q and U in the 
Auriga region are somewhat steeper than in Horologium, 
indicating that the Horologium region probably contains 
more small-scale structure in the Faraday screen than the 
Auriga region. 

The power spectra of RM in Fig.0are shallower than 
the Q, U or P power spectra. In fact, we do not expect 
a direct correspondence between the multipole spectral 
indices of RM and P (or Q, U), as the former describes 
very directly the electron content and magnetic field in 
the ISM (integrated over the line of sight), whereas in 
the latter case the polarized radiation is modulated by 
Faraday rotation and depolarization. 

3.3. Power spectra from the WENSS polarization 
region 

In the WENSS polarization region, power spectra were 
evaluated for subfields, to study possible dependences of 
the multipole spectral index on Galactic longitude and/or 
latitude. The 11 subfields are shown in Fig. ED superim- 
posed on grey scale maps of P. The power spectra of polar- 
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Table 1. Multipole spectral indices a for observed polarized intensity P, Stokes Q and Stokes U, for 5 frequencies 
and their average over frequency, and olrm, m the Auriga and Horologium regions. Values for otp, oiq and a\j are 
given for tapered data, denoted by a subscript H\ and untapered data. The multipole ranges used to derive a were 
400 < I < 1500. Only for ap of the untapered data, the range was smaller because of flattening of the spectrum at 
higher I. 



Auriga 





341 MHz 


349 MHz 


355 MHz 


360 MHz 


375 MHz 


mean 


ap,t 


2.37 ±0.19 


2.59 ±0.20 


2.38 ±0.19 


2.35 ±0.17 


1.88 ±0.18 


2.32 ±0.08 


ap 


2.06 ±0.54 


2.23 ±0.46 


2.49 ±0.60 


2.00 ±0.47 


2.20 ±0.59 


2.20 ±0.24 


a Q,t 


3.55 ±0.20 


3.72 ±0.21 


3.59 ±0.20 


3.43 ±0.19 


2.50 ±0.20 


3.36 ±0.09 


olq 


3.13 ±0.24 


3.16±0.25 


3.16 ±0.24 


3.00 ±0.24 


3.12±0.25 


3.12±0.11 


uu,t 


3.69 ±0.19 


3.89 ±0.21 


3.79 ±0.21 


3.69 ±0.21 


3.49 ±0.20 


3.71 ±0.09 


OLJJ 


3.25 ±0.24 


3.20 ±0.25 


3.23 ±0.24 


3.13 ±0.24 


3.18 ±0.53 


3.20 ±0.11 


OLRM,t 












0.99 ±0.08 








Horologium 








341 MHz 


349 MHz 


355 MHz 


360 MHz 


375 MHz 


mean 


ap,t 


2.11±0.19 


2.18±0.19 


2.11±0.18 


1.98 ±0.19 


1.95 ±0.19 


2.07 ±0.08 


ap 


2.12 ±0.63 


2.36 ±0.55 


2.26 ±0.66 


2.47 ±0.68 


2.49 ±0.91 


2.34 ±0.31 


a Q,t 


2.58 ±0.22 


3.17±0.19 


2.96 ±0.19 


2.90 ±0.20 


2.78 ±0.20 


2.89 ±0.09 


olq 


2.15 ±0.26 


2.92 ±0.22 


2.57 ±0.23 


2.60 ±0.23 


2.32 ±0.25 


2.51±0.11 


uu,t 


2.56 ±0.22 


3.14 ±0.20 


3.04 ±0.19 


3.22 ±0.20 


3.24 ±0.17 


3.04 ±0.09 


a v 


2.07 ±0.26 


3.00 ±0.22 


2.69 ±0.23 


2.85 ±0.23 


2.65 ±0.25 


2.65 ±0.11 


CXRM,t 












0.94 ±0.10 



w 
a, 
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Fig. 4. Power spectrum of P at 341 MHz in the Auriga 
field. The solid line shows the observed data, the dotted 
line is corrected for the 300m taper. 
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ized intensity P are shown in Fig. 1101 where the subfields 
are arranged as in Fig. |5J The power spectra in subfields 
9, 10 and 11, at high Galactic latitude b, have a lower 
amplitude than the power spectra at lower b, which is 
consistent with the decreasing amount of P at higher b 
visible in Fig. 

The multipole spectral indices of the power spectra of 
P, Q and U are given in Table and the dependence of 
ap on Galactic longitude and latitude is shown in Fig. ^2 



Fig. 5. Power spectra of polarized intensity P for 5 fre- 
quency bands in the Auriga region (top) and Horologium 
region (bottom). In each plot, the upper line of symbols 
denotes the tapered data, the lower line the untapered 
data, and the solid lines are linear fits to the spectra. In 
the untapered data, the spectrum is flattened at high I 
due to noise. 

The observed decrease of spectral index with increasing 
latitude (i.e. power spectra become flatter with increas- 
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Fig. 6. Power spectra of RM in the Auriga region (top) 
and Horologium region (bottom), derived from the ta- 
pered data. 
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Fig. 7. Power spectra of Stokes Q (top) and Stokes U 
(bottom) for 5 frequencies in the Auriga region. Notation 
as in Fig. [5] 



ing latitude) indicates a decrease in the amount of large- 
scale structure with increasing latitude. The dependence 
of spectral index on longitude is computed after rescaling 
of the data to a standard latitude of b — 15°. The fitted 
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Fig. 8. Power spectra of Stokes Q (top) and Stokes U (bot- 
tom) for 5 frequencies in the Horologium region. Notation 
as in Fig. 03 

Table 2. Multipole spectral indices a for observed polar- 
ized intensity P, Stokes Q and Stokes U in 11 subfields in 
the WENSS polarization region at 325 MHz. 



No. 


(l,b) (V) 


ap 


OiQ 


ecu 


1 


(159, 04) 


1.67 ±0.08 


2.55 ±0.09 


2.43 ±0.10 


2 


(165, 11) 


1.48 ±0.08 


2.10 ±0.08 


2.06 ±0.08 


3 


(158, 11) 


1.84 ±0.09 


2.50 ±0.09 


2.20 ±0.09 


4 


(151, 11) 


1.70 ±0.08 


2.09 ±0.09 


2.02 ±0.09 


5 


(144, 11) 


1.63 ±0.08 


2.28 ±0.10 


2.31 ±0.10 


6 


(143, 18) 


1.24 ±0.09 


1.96 ±0.10 


2.04 ±0.10 


7 


(150, 18) 


0.99 ±0.08 


1.44 ±0.08 


1.47 ±0.09 


8 


(157, 18) 


1.67 ±0.09 


2.33 ±0.09 


2.33 ±0.09 


9 


(157, 25) 


0.73 ±0.08 


1.24 ±0.08 


1.18 ±0.09 


10 


(150, 25) 


0.90 ±0.08 


1.18 ±0.08 


1.44 ±0.08 


11 


(144, 25) 


0.68 ±0.08 


0.77 ±0.08 


0.63 ±0.08 



slopes in Fig. 1111 and slopes determined for cxq and ajj in 
a similar way, are: 

dap/db = -0.051 ± 0.004 da P /dl = -0.001 ± 0.003 
da Q fdb = -0.073 ± 0.004 da Q jdl = -0.008 ± 0.004 
dau/db = -0.066 ± 0.004 da v /dl = -0.005 ± 0.004 

In summary, the spectral index a in P, Q and U decreases 
with increasing Galactic latitude, and is consistent with no 
dependence of a on Galactic longitude. Although the small 
errors in the derived slopes suggest a good determination 
of the slope, the large spread of the data points in Fig. 1111 
indicates that a linear gradient is not the perfect model 
to describe the data. 

The decrease of a with latitude indicates that there is 
more structure in polarization on larger scales at lower lat- 
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Fig. 11. Dependence of multipole spectral index ap on 
Galactic latitude b (left) and longitude I (right), for 11 
subfields in the polarization part of the WENSS survey. 
The spectral indices in the right plot are normalized to 
the latitude gradient shown in the left plot. 




i n 

Fig. 9. Polarized intensity at 325 MHz of the WENSS 
region in grey scale, superimposed with white boxes de- 
noting the numbered subfields over which angular power 
spectra were computed. 
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Fig. 10. Power spectra of polarized intensity P in sub- 
fields in the WENSS polarization region. The plots are 
arranged as in Fig. 



itude. Thus, although small-scale structure in polarization 
is seen up to very high latitudes, even at frequencies as low 
as 350 MHz (Katgert & de Bruyn 1999), the amount of 
structure and its spectral index decrease with latitude, at 
least in this region of the sky. The region studied here may 
be special, as it is situated at the edge of a region of high 
polarization (the "fan region" , see e.g. Brouw & Spoelstra 
1976), which is thought to have little structure in magnetic 
field and/or electron density compared to its surroundings. 
This high-polarization region extends to b « 20°, and its 
edge could be responsible for the decrease in structure on 
scales of ~ 1°. 



The spectral indices in the Auriga (b = 16°) and 
Horologium (b = 7°) regions are higher than implied by 
Fig. 1111 This might be due to the way in which the Auriga 
and Horologium regions were selected, viz. because of their 
remarkable structure on degree scales. 

3.4. Existing literature of power spectra from diffuse 
polarization 

Much work has been done on the determination of power 
spectra of the diffuse Galactic synchrotron background, 
because the Galactic synchrotron radiation is a fore- 
ground contaminator in Cosmic Microwave Background 
Radiation (CMBR) polarization measurements at high 
frequencies v ss 30 - 100 GHz. Power spectra of the diffuse 
polarized synchrotron background intensity have been de- 
termined from several radio surveys at frequencies from 
408 MHz to 2.7 GHz, in many parts of the sky (Tucci 
et al. 2000, 2002, Baccigalupi et al. 2001, Giardino et al. 
2002, Bruscoli et al. 2002). 

These power spectra studies are based on the following 
surveys of polarized radiation: 

— Dwingeloo 25m-dish survey (Brouw & Spoelstra 1976), 
of the region 120° < I < 180° and b > -10°(although 
undersampled). This is a multi-frequency survey at 
408 MHz, 465 MHz, 610 MHz, 820 MHz and 1411 MHz, 
with increasing angular resolution of 2.3° to 0.5°. 

— Parkes 2.4 GHz Galactic plane survey (Duncan et al. 
1997), of the region 238° < I < 5° and with \b\ < 5°, 
at some positions a few degrees higher, at a resolution 
of 10.4'. 

— Effclsberg 2.695 GHz Galactic plane survey (Duncan 
et al. 1999), of the region 5° < I < 74°, and |6| < 5°, 
at a resolution of 4.3'. 

— Effelsberg 1.4 GHz intermediate latitude survey 
(Uyamker et al. 1999), which consists of 4 regions 
within 45° < I < 210° and -15° < b < 20° with 
an angular resolution of 10.4'. 

— Australia Telescope Compact Array (ATCA) 1.4 GHz 
survey (Gaensler et al. 2001). This is a test region for 
the Southern Galactic Plane Survey (SGPS, McClure- 
Griffiths et al. 2001) at 325.5° < I < 332.5°, -0.5° < 
b < 3.5°, with a resolution of about 1'. 
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Power spectra of total intensity / and polarized intensity 
P were derived in these surveys for multipoles over a range 
of i « 10 to 6000. 

Fig. 1121 shows the variation of ap with Galactic longi- 
tude, latitude and frequency, using the available data as 
detailed in Table |3| In the left plots, the lines show ranges 
in longitude (top) and latitude (bottom) over which ap 
was computed. Solid lines give high multipolc numbers 
(100 < £ < 6000), dashed-dotted lines denote an inter- 
mediate multipole range (30 < £ < 200) and the dotted 
lines give small multipoles (10 < t < 80). The WSRT data 
from the Auriga, Horologium en WENSS fields, discussed 
here, are given in asterisks. In the right plot, ap against 
frequency is displayed. The connected points are from the 
same area observed at different frequencies from 408 MHz 
to 1411 MHz (in the Brouw & Spoelstra paper). 

The first conclusion from Fig.EJis that multipole spec- 
tral indices vary over the sky from a sa 1 to 3, without 
showing a clear correlation with Galactic longitude, while 
only the WENSS subfields show a dependence of spectral 
index on latitude. However, all surveys have been done at 
different resolutions and frequencies, and the regions used 
to compute power spectra are of different sizes. This can 
explain why a possible dependence of a on latitude was 
not clearly seen in the other studies. The large variation 
in slopes of angular power spectra in P indicates that in- 
terpretation of the slope is not straightforward, possibly 
due to large influence of depolarization mechanisms. Care 
must therefore be taken in extrapolating the results to 
higher frequencies. 

Furthermore, spectral indices show an increase with 
frequency from 408 MHz to 1.4 GHz. This means that 
the power spectra become steeper, so the relative amount 
of small-scale structure decreases. This could be due to 
the large Faraday rotation at low frequencies. Typical 
RMs of 5 rad m~ 2 are present in the Brouw & Spoelstra 
data (Spoelstra 1984), and will rotate polarization an- 
gles at 325 MHz by about 250°. Variations in RM of 
a few rad m -2 give angle variations of over 90°, which 
would cause beam depolarization if the angle variations 
occur on scales smaller than the beam (2.3° in this case). 
Beam depolarization only acts on scales of the synthesized 
beam, and therefore creates additional structure on small 
scales in P, which flattens the power spectrum. A ARM 
of 5 rad m~ 2 would cause a variation in polarization an- 
gle of about 40° at 820 MHz, and of no more than 10° 
at 1.4 GHz. So at frequencies above 1.4 GHz, a ARM of 
5 rad m -2 would cause negligible beam depolarization. In 
addition, the resolution of the observations generally in- 
creases with increasing frequency, which would also cause 
a decrease in beam depolarization. This might explain why 
above 1.4 GHz the spectral index does not appear to be 
correlated with frequency. The fact that spectral indices 
in the two WSRT regions are much higher than would be 
expected from this argument can be due to the criteria 
used to select the two fields. 



4. Structure functions 

The disadvantage of using angular power spectra is that 
a regular grid of data is required. If the data are very ir- 
regularly spaced (e.g. in the case of data from pulsars or 
extragalactic point sources), it is better to use the struc- 
ture function which in principle gives the same informa- 
tion, but can be calculated easily for irregularly spaced 
data. The structure function SF of a radiation field X, as 
a function of distance lag d, is 

SFx(d) _ skw*) (2) 

where X(xi) is the value of field X at position xi and N is 
the number of data points. If the power spectrum of RM 
is a power law with spectral index a, then the structure 
function SFrm is 

on i j\ iu i f a - 2 for 2 < a < 4 

SF RM (d) oc eP where M = S 2 for a > 4 ^ 

(Simonetti et al. 1984). We determined the structure func- 
tions of RM to compare with existing estimates of the 
structure function of Galactic RM from polarized extra- 
galactic point sources. As the determination of structure 
functions does not require a regular grid, we can compute 
the SF including and excluding "bad" data points to ex- 
amine how the structure functions of P and RM change. 
This will allow us to estimate the effect of "bad" data in 
the power spectra in P and RM computed in Sect.|3| 

4.1. Structure functions of RM 

Structure functions SFrm in the Auriga and Horologium 
regions are plotted against distance d in degrees in the log- 
log plots in Fig. 1131 (solid line), where error bars denote 
the standard deviation. The minimum distance shown is 
d w 7'. For the evaluation of the SF, only "reliably deter- 
mined" RM values are used, according to the definition 
in Sect. 12.11 Although the spectrum is consistent with a 
flat slope, there is some evidence for a break in the slope 
at d — 0.3°, primarily in the Horologium field. For larger 
angular scales, the SF is approximately flat in the Auriga 
region, with a tentative increase at the largest lags, and 
even decreasing in Horologium. 

We can estimate the magnitude of the contribution of 
unreliably determined RMs by reevaluating the SF for the 
complete grid of RM values, instead of only the reliable 
RMs. This estimate is important for a discussion of the 
power spectra of RM , which were evaluated over the com- 
plete dataset, including unreliably determined RMs. The 
structure function using the complete dataset is shown 
in the left panel of Fig. as a dotted line. The struc- 
ture function clearly has a lower amplitude if the unreli- 
able RM determinations are removed, but the slope of the 
structure function remains approximately the same. 
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Fig. 12. Estimates of multipole spectral indices ap(£) from our observations and from the literature, as a function of 
Galactic longitude (top), Galactic latitude (bottom), and frequency (right). In the left-hand plots, solid lines denote 
high multipole numbers (100 < I < 6000), the dashed-dotted line is an intermediate multipole range (30 < I < 200) 
and the dotted lines give small multipoles (10 < I < 80). In the right-hand (frequency) plot, the dotted lines connect 
observations of the same region made at different frequencies. The WSRT observations discussed in this paper are 
shown by asterisks. 
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Fig. 13. The solid lines show the RM structure function 
SFrm as a function of distance d in degrees for the Auriga 
region (top) and the Horologium region (bottom), using 
only reliable RM values. The dotted line in the top plot 
gives SFrm evaluated for the entire grid, including unre- 
liable RM determinations. 



4.2. Structure functions of P 

We compute the structure function of polarized inten- 
sity, for both the complete grid of beams, and for those 
beams selected to have high P. In Fig. we show struc- 
ture functions of P in the Auriga region for 5 frequencies, 
where again the error bars denote the standard deviation. 
The average logarithmic slope of the structure functions 
is ~ 0.35 in the range 0.2° ^ d 1°, and the spectrum 
flattens on larger scales, probably due to the insensitivity 
of the WSRT to structure on large angular scales. The 
solid line denotes the estimate based only on beams with 
P > 4a, and the dotted line gives the estimate based on all 
the data. The structure functions based on selected beams 
and those based on all data are not significantly different 
on small scales, and start deviating only for d J; 0.7°. 
Although small deviations are created by selecting the 
best data, the overall slope of the structure function is 
hardly affected by the selection including only regions of 
high P. The relatively small effect of including bad data 
gives confidence that the determination of the logarithmic 
slope of the power spectra of RM and P is not significantly 
influenced by noisy data or poorly determined RM. 

5. Discussion 

To the authors knowledge, the results presented here are 
the first observational determinations of the power spec- 
trum of rotation measure. What power spectrum is to be 
expected is unclear, as the RM is a complex quantity de- 
pending on magnetic field structure and direction, elec- 
tron density and the path length through the ionized ISM. 
There are indications that the electron density and RM 
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Fig. 14. Structure function SFp as a function of distance d in degrees for the Auriga region, in 5 frequencies. The 
solid line shows SFp in which only beams with P > 4<r are included, the dotted line gives SFp as computed from all 
data. 



exhibit Kolmogorov turbulence (Minter & Spangler 1996, 
Armstrong et al. 1995). 

An approximately flat structure function corresponds 
to structure of equal amplitude on all scales, similar to 
a noise spectrum. The break in SFrm at 89 w 0.3°, if 
present, corresponds to a change in characteristics of the 
structure on scales of ~ 3.9 pc, assuming a path length of 
600 pc (Haverkorn et al. 2003d). 

Minter & Spangler (1996) studied structure functions 
of RM from polarized extragalactic sources and of emis- 
sion measure (EM) from Ha measurements, on the same 
angular scales as in our observations. They find a break 
in the slope of the structure function, which can be in- 
terpreted as a transition of 3D Kolmogorov turbulence 
(fi = 5/3, where fi is the slope of the structure function 
as defined in Eq. ©) to 2D turbulence (/i = 2/3) in RM 
on angular scales of 89 w 0.1°, the scale of our resolution. 
They assume a total path length of 2000 pc, and conclude 
that the transition occurs at scales of ~ 3.6 pc, consistent 
with the spatial scale at which we estimate the turnover 
point in the SF of the Auriga and Horologium regions. 
However, there is a large discrepancy between fi values 
found by Minter & Spangler (/z = 5/3 and 2/3) and by us 
(fi w 0), which could be explained by the fact that they 
probe the complete line of sight through the medium of 
many kpc, whereas the RM that we obtain is only pro- 
duced in the nearest ~ 600 pc. Therefore, it is possible 
that the nature of the turbulence changes from /i = 
nearby (approximately in the Galactic stellar disk), to 2D 
and/or Kolmogorov-like turbulence at larger distances. 

However, flat structure functions are also found by 
Simonetti et al. (1984), who have studied structure func- 
tions from RMs in three regions of the sky. Two of the 
three regions (around the North Galactic Pole and at 
180° < I < 220°, 10° < b < 50°) are consistent with a flat 
structure function (/z = 0). (The third region is located 
along the Local arm and in the plane, at 70° < I < 110° 
and —45° < b < 5, and its structure function of RM has 
a slope /j w 1.) The smallest angular scale they probe in 



these regions is about 2°. Clegg et al. (1992) study low- 
latitude extragalactic sources and also obtain flat struc- 
ture functions, albeit with a much higher amplitude than 
that of the high-latitude extragalactic sources. 

6. Conclusions 

The multipole spectral index for polarized intensity P is 
a ~ 2.2 in the Auriga and Horologium regions, and ranges 
from 0.7 to 1.8 for subfields in the WENSS polarization re- 
gion. The multipole spectral index decreases with Galactic 
latitude (i.e. power spectra become flatter towards higher 
latitudes), but is probably constant with Galactic longi- 
tude. 

In all regions, the power spectra of Stokes Q and U 
are steeper than the spectra of P. This is most likely due 
to the presence of a Faraday screen, which creates addi- 
tional structure in Q and U, but not in P. As the Faraday 
screen is located in front of the polarized emission, the 
structure induced by the screen will be on larger angu- 
lar scales than that of the emission, which steepens the 
Q and U spectra. The derived power spectra of P agree 
with earlier estimates, although all estimates show a large 
range of 0.6 ap Sa 3, possibly due to a large influence 
of depolarization mechanisms. This makes interpretation 
of the power spectra uncertain. 

Structure functions of RM in the Auriga and 
Horologium fields are consistent with flat spectra (i.e. a 
logarithmic slope of the structure function /Lt = 0), but 
may show a break close to 0.3°, which is at the same spa- 
tial scales as a break in the structure function in the RMs 
of extragalactic sources (Minter & Spangler 1996). The 
flat spectrum indicates a noise-like spectrum with equal 
amounts of structure on all scales. 

The derived structure functions support the estimates 
of the power spectra, because they are based only on high- 
quality data, i.e. reliably determined RM, and P with 
high S/N. Structure functions using all data do not show 
a significant difference with the structure functions of the 
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selected data. This gives confidence that the power spectra 
determinations are not significantly affected by bad data 
points. 
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280 - 


290 


11 




11 


" 




1.91 






290 - 


300 


11 




11 


" 




1.77 






300 - 


310 


11 




11 


" 




1.28 






310 - 


320 


11 




11 






1.64 






320 - 


330 


11 




11 






1.93 






330 - 


340 


11 




11 


" 




1.44 






340 - 


350 


11 




11 


" 




1.47 






350 - 


360 


11 




11 


" 




1.49 






240 - 


360 


11 




11 


" 




1.7 


Bruscoli et al. 2002 




326.5 - 


331.5 


- 1 - 


4 


11 


" 




1.68 


Tucci et al. 2002 




240 - 


360 


- 5 - 


5 


11 


40 - 


250 


2.4 


Giardino et al. 2002 


Effelsberg 


20 - 


30 


- 5 - 


5 


2695 


100 - 


800 


1.79 


Baccigalupi et al. 2001 


Duncan et al. 


30 - 


40 


15 




51 


51 




1.72 




(1999) 


40 - 


50 


11 




11 


11 




1.55 






50 - 


60 


11 




11 


11 




1.98 






55 - 


65 


11 




11 


11 




1.93 






5 - 


75 


11 




11 


11 




1.6 


Bruscoli et al. 2002 


Effelsberg 


45 - 


55 


5 - 


20 


1400 


100 - 


800 


1.5 


Bruscoli et al. 2002 


Uyaniker et al. 


140 - 


150 


4 - 


10 


51 


51 




2.5 




(1999) 


190 - 


200 


8 - 


15 


51 


51 




2.3 





Table 3. Studies of power spectra of diffuse polarized intensity from different surveys. Given are the range in Galactic 
longitude I and latitude b, the frequency of observation v, the range in multipolc I for which the power spectrum was 
computed, the spectral index of P ap and the reference to the study of the power spectra. 



M. Haverkorn et al.: Structure in the Galactic polarized radio background 



survey 




) 


b(° 


) 


v (MHz) 




t 


Ctp 


reference 


ATCA 


327 - 


331 


- 0.5 - 


3.5 


1400 


600 - 


6000 


1.68 


Tucci et al. 2002 


Gaensler et al. 


329 - 


332 


- 


3 


» 






1.66 




(2001) 


326 - 


329 


11 










1.68 






330.1 - 


331.1 


0.5 - 


1.5 


» 






1.79 






329.2 - 


330.2 


0.6 - 


1.6 


" 






2.24 






327.2 - 


328.2 






" 






1.93 






328 - 


329 


2.1 - 


3.1 


" 






2.58 






326.5 - 


327.5 


1.8 - 


2.8 


" 






1.88 






329.9 - 


330.9 


It 




" 






1.68 




WSRT 


158 - 


165 


13 - 


20 


350 


100 - 


1000 


2.26 


Haverkorn et al. 2003a 


WSRT 


134 - 


141 


3 - 


10 


350 


100 - 


1000 


2.20 


Haverkorn et al. 2003b 


WSRT 


156 - 


143 


0.5 - 


7.5 


327 


100 - 


1500 


1.67 


Schnitzeler et al. in prep 




162 - 


169 


8 - 


15 


» 






1.48 






155.5 - 


161.5 


" 










1.84 






147.5 - 


154.5 


" 










1.70 






140.5 - 


147.5 


» 




" 






1.63 






139 - 


146 


15 - 


22 


" 






1.24 






146.5 - 


153.5 


" 










0.99 






153.5 - 


160.5 












1.67 






153.3 - 


160.5 


22 - 


29 








0.73 






147 - 


154 












0.90 






140.5 - 


147.5 












0.68 





Table 3. Continued 



